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Abstract 


A  theory  is  developed  for  the  incorporation  of  conduc¬ 
tion  band  anisotropy  into  the  analysis  of  the  excitcn-plasma 
Mott  transition  in  indirect  gap  semiconductors.  Ellipsoidal 
energy  surfaces  are  assumed  for  the  electrons  while  spheri¬ 
cal  energy  surfaces  are  retained  for  holes.  Static  electron 
hole  screening  in  the  random  phase  approximation  is  assumed. 
The  Mott  transition  is  associated  with  the  electron-hole 
pair  density  at  which  the  exciton  binding  energy  in  the 
assumed  potential  is  zero.  The  binding  energy  is  computed 
^ariationedly . 

It  is  found  that  the  electron  anisotropy  causes  the 
MOtt  transition  to  shift  to  higher  densities.  It  is  also 
^ound  that,  in  the  absence  of  screening,  the  exciton  binding 
energy  is  not  significantly  affected  by  the  electron  aniso¬ 
tropy.  It  is  thus  concluded  that  the  shift  to  higher  densi¬ 
ties  is  due  Icurgely  to  the  reduced  ability  of  emisotropic 
electrons  to  screen. 


EFFECTS  OF  CONDUCTION  BAND  ANISOTROPY  ON  THE  EXCITON-PLASMA 
MOTT  TRANSITION  IN  INDIRECT  GAP  SEMICONDUCTORS 

% 

I.  Introduction 

The  exciton-plasma  Mott  transition  in  silicon  has  been 
studied  in  detail  in  receht  years.  The  experimental  data 
(Forchel,  1982;  Forchel,  to  be  published)  has  been  success¬ 
fully  explained  (Norris  and  Bajaj,  1982)  by  assuming  that 
the  electron-hole  interaction  is  statically  screened,  and  by 
taking  the  conduction  and  valence  bands  to  be  isotropic. 

While  the  assximption  of  conduction  band  isotropy  works  well 
in  Si,  where  the  ratio  of  longitudinal  to  transverse  electron 
masses  is  less  than  five,  one  would  expect  the  effects  of 
conduction  band  amisotropy  to  be  larger  in  Ge,  where  the  ratio 
of  longitudinal  to  transverse  electron  masses  is  greater  than 
nineteen.  It  is  thus  the  purpose  of  this  work  to  extend  the 
previously  mentioned  theory  by  taking  conduction  band  ani¬ 
sotropy  into  account. 

The  system  under  study  here  consists  of  optically 
generated  electron-hole  (E-H)  pairs  in  pure  indirect  gap 
semiconductors.  The  presentation  will  therefore  begin  with 
a  discussion  of  how  excitons  are  formed  through  optical 
pumping,  followed  by  the  definition  of  the  exciton-plasma 
Mott  transition. 

In  the  next  chapter,  the  experimental  work  will  be 
discussed.  This  will  be  done  to  show  the  experimental 


evidence  for  the  occurrence  of  the  exciton-plasma  Mott 
transition,  and  to  show  clearly  how  the  connection  between 
theory  and  experiment  is  made. 

The  theory  will  be  developed  in  chapter  three  and 
results  for  both  Si  and  Ge  will  be  presented  in  chapter 
four.  Finally,  the  conclusions  will  be  summarized  in 
chapter  five. 

Exciton  Formation 

The  formation  of  excitons  in  an  indirect  gap  semicon¬ 
ductor  is  shown  schematically  in  Figure  1  (Wolfe,  1982) .  A 
photon,  which  has  an  energy  greater  th2ui  the  energy  gap  (Eg) 
of  the  semiconductor,  is  absorbed  emd  excites  an  electron 
from  the  valence  band  to  a  conduction  band  state  which  lies 
above  the  conduction  band  minimum.  The  electron  then  under¬ 
goes  a  rapid  thermalization  process  in  which  it  loses  energy 
to  the  lattice  through  the  emission  of  phonons  euid  thus 
relaxes  to  the  conduction  bemd  minimum.  Hydrogen-like 
bound  states  exist  within  the  energy  gap.  They  are  due  to 
bound  E-H  pairs,  which  are  called  excitons,  zmd  which  may 
form  if  the  sample  is  sufficiently  cool  (e.g.  T  <  80*K  for 
Si) .  The  electron  and  hole  eventually  recombine  giving  off 
a  characteristic  luminescence.  Since  the  material  is  in¬ 
direct  gap,  the  E-H  recombination  is  also  accompanied  by 
the  emission  or  absorption  of  phonons. 


Exclton  Formation  (Wclfe,  1982) 


The  Exciton-Plasma  Mott  Transition 


Consider  a  free  exciton  (FE)  gas  and  ask  the  question, 
"How  can  the  bound  electrons  and  holes  become  dissociated?" 

One  way  in  which  the  E-H  pairs  may  become  dissociated 
is  through  thermal  ionization.  This  is  a  diffuse  process  in 
that,  for  a  given  temperature,  the  FE  gas  and  EHP  coexist  in 
thermal  equilibrium  and  there  is  not  a  precisely  defined  E-H 
pair  density  at  which  dissociation  occurs. 

Another  way  in  which  E-H  pairs  may  become  dissociated 
is  through  entropy  ionization  (see  for  example;  Mock,  1978) . 
This  is  a  complicated  effect  which  occurs  when  the  total 
number  of  electrons  and  holes  in  the  system  is  reduced  at. 
constant  temperature.  The  important  point  here  is  again 
that  there  is  no  precisely  defined  E-H  pair  density  at  which 
dissociation  occurs. 

The  third  way  in  which  excitons  may  ionize  is  through 
screening.  In  this  case  the  exciton  density  is  increased, 
at  constant  temperature,  by  increasing  the  optical  pumping 
power  and  a  point  is  reached  at  which  a  given  electron  can 
no  longer  be  associated  with  any  particular  hole.  Here  there 
is  a  precisely  defined  density  at  v/hich  the  exitonic  binding 
energy  goes  to  zero.  The  transition  is  from  an  insulating 
FE  gas  to  a  metallic  electron-hole  plasma  (EHP)  and  is  called 
a  Mott  transition. 

Thus,  the  Mott  transition  is  an  insulator-to-metal 
transition  which  arises  as  a  result  of  screening. 
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II .  Experiment 


Experimentally,  the  exci ton-plasma  Mott  transition  has 
been  studied  in  connection  with  a  separate  transition  which 
is  known  to  occur  and  which  is  well  understood.  This  other 
transition  is  a  first-order  phase  transition  in  which  the 
PE  gas  condenses  into  a  highly  dense  metallic  electron-hole 
liquid  (EHL) .  Experimental  work  may  be  seen  for  example  in 
Hammond  et  al  (Hammond,  1976) ,  Thomas  et  al  (Thomas,  1973) , 
emd  Thomas  et  al  (Thomas,  1974).  A  thorough  discussion  of 
the  theory  of  the  EHL  is  given  by  Rice  (Rice,  1977) . 

The  EHL  is  important  in  the  experimental  study  of  the 
Mott  transition  and  will  thus  be  briefly  discussed. 

The  phase  diagram  for  the  electron-hole  system  is  shown 
in  Figure  2.  How  such  a  phase  diagram  is  constructed  will 
be  discussed  shortly.  For  the  present  it  is  sufficient  to 
note  that  below  a  certain  critical  temperature,  T^,,  a  two- 
phase  region  exists  where  the  FE  gas  (or  EHP)  is  in  equili- 
britun  with  the  dense  EHL.  The  solid  curve  is  the  liquid-gas 
coexistence  curve.  If  the  FE  gas  density  is  increased  (by 
increased  optical  pumping  at  sufficiently  low,  constant 
temperature) ,  the  gas-to-liquid  transition  will  occur  when 
the  FE  gas  density  reaches  the  value  on  the  coexistence 
curve.  The  EHL  will  then  be  present  in  droplets  whose  den¬ 
sity  is  determined  by  the  liquid  side  of  the  coexistence 
curve.  Further  increases  in  optical  pumping  change  the  size 
of  the  electron-hole  droplets  but  do  not  change  their  density 


Given  the  above  comments  on  the  EHL,  one  may  now  turn 
to  the  experimental  construction  of  the  phase  diagram  of 
Figure  2.  This  construction  will  be  illustrated  through  a 
discussion  of  a  specific  experimental  study  of  the  exciton- 
plasma  Mott  transition  (Shah,  1977) .  Shah  emd  his  co-workers 
were  first  to  report  experimental  evidence  for  the  Mott  dis¬ 
sociation  of  excitons  into  EHP.  Their  work  shows  clearly 
how  the  experimental  data  are  obtained  and  interpreted.  This 
discussion  will  thus  serve  to  define  the  theoretical  problem 
as  well  as  to  provide  the  experimental  background. 

The  Luminescence  Spectra 

Shah  and  his  co-workers  excited  a  crystal  of  pure  Si 
with  radiation  from  an  argon  laser  and  observed  the  resulting 
luminescence  spectra  at  various  temperatures  and  for  various 
optical  pumping  powers. 

The  spectra  which  they  observed  at  low  temperature  (18*K) 
are  shown  in  Figure  3.  For  low  pumping  powers,  the  FE  line 
was  seen,  while  for  sufficiently  high  pumping  powers,  a  second 
peak  due  to  EHL  Ivuninescence  was  also  present.  Both  of  these 
luminescence  peaks  were  found  to  have  line  shapes  which  were 
independent  of  pumping  power. 

Figure  4  shows  the  spectra  which  were  obtained  at  23 *K. 
For  low  pumping  power,  the  FE  line  was  again  seen.  However, 
as  the  incident  power  was  increased,  the  low  energy  side  of 
the  FE  line  was  observed  to  broaden  as  shown  by  the  singly 
dashed  curve.  Again,  for  sufficiently  high  pumping  powers. 
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18’K 


1.12 
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the  second  peak  due  to  the  formation  of  EHL  v/as  observed. 


The  spectra  obtained  at  30*K,  which  is  above  the  criti¬ 
cal  temperature  for  the  formation  of  an  EHL,  are  shown  in 
Figure  5.  At  low  incident  power  the  FE  line  was  observed 
(curve  1) ,  and  as  the  power  was  increased,  the  line  broaden¬ 
ing  was  seen  to  evolve  continuously  into  a  shape  which  is 
well  fit  by  a  plasma  line  shape  (curve  3) . 

Shah  emd  his  co-workers  interpreted  the  observed  line 
broadening  as  being  due  to  the  Mott  dissociation  of  excitons 
into  an  EHP.  This  is  a  logical  interpretation  since  the  FE 
line  shape  is  independent  of  pumping  power,  while  the  line 
shape  due  to  the  recombination  of  unboiind  electrons  and  holes 
does  depend  on  pumping  power.  It  is  thus  natural  to  associate 
the  onset  of  the  line  broadening  with  the  onset  of  the  Mott 
dissociation.  How  this  was  quantified  will  now  be  presented. 

The  Threshold  Pumping  Powers 

In  order  to  determine  the  threshold  powers  for  the  Mott 
transition,  and  for  the  formation  of  EHL,  Shah  and  his  co¬ 
workers  plotted  the  change  in  position  of  the  low  energy 
half-maximum  of  each  distinct  peak  observed  in  the  spectra. 

The  ch2uige  was  plotted  as  a  function  of  average  incident 
power.  The  results  are  shown  in  Figure  6. 

At  18 '*K,  only  the  FE  and  EHL  peaks  were  seen  and  their 
shapes  were  independent  of  pumping  power.  Thus  the  power, 

1^,  at  which  the  EHL  peak  was  first  observed  was  taken  as 
the  threshold  for  EHL  formation. 


{.30*K)  (Shah,  1977) 


At  25 *K,  the  line  broadening  was  found  to  yield  a  linear 
plot  between  the  FE  and  EHL  plots.  was  determined  as 
before  for  EHL  condensation.  The  threshold,  1^,  for  the 
onset  of  line  broadening  was  determined  by  extrapolation  as 
shown. 

At  30 the  continhous  broadening  was  observed  and  Iq 
was  again  determined  by  extrapolation. 

Once  the  threshold  powers  were  determined,  there 
remained  the  problem  of  converting  those  powers  into  E-H 
pair  densities.  This  problem  was  solved  by  constructing  the 
phase  diagram  for  the  E-H  system. 

The  Phase  Diagram 

The  phase  diagram  for  the  E-H  system  in  Si  was  shown  in 
Figure  2.  The  solid  curve  is  taken  from  Norris  and  Bajaj 
(Norris,  1982) ,  while  the  circle  and  triauigle  points  are 
Included  to  aid  in  the  description  of  how  Shah  et  al  con¬ 
structed  a  similar  phase  diagram. 

The  EHL  densities  (triangular  points)  were  determined 
from  the  liquid  luminescence  half-width  by  using  the  theo¬ 
retical  calculations  of  Hammond,  McGill,  and  Mayer  (Heunmond, 
1976) .  This  process  resulted  in  the  experimental  determina¬ 
tion  of  n^,  the  EHL  density  at  0®K. 

Given  experimental  values  for  n^^  and  T^  (the  critical 
temperature  for  EHL  condensation)  the  theoretical  calcula¬ 
tions  of  Reinecke  and  Ying  (Reinecke,  1975)  were  used  to 
obtain  the  liquid-gas  coexistence  curve  (solid  line) .  It 
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was  thus  possible  to  determine  the  FE  gas  densities  (circles 
on  the  solid  curve)  which  were  in  equilibrium  with  the  EHL 
at  various  temperatures,  and  to  associate  these  densities 


with  the  measured  threshold  powers  (I^)  for  EHL  condensation 

The  Mott  transition  densities  were  then  determined  from 
the  threshold  powers,  1^,’  by  assviming  a  temperature-indepen¬ 
dent,  linear  scaling  between  E-H  pair  density  and  pvunping 
power.  This  assumption  is  of  uncertain  validity  but  repre¬ 
sents  the  best  method  available  for  the  experimental  deter¬ 
mination  of  the  Mott  transition  densities. 

The  Theoretical  Problem 

The  above  discussion  clearly  defines  the  theoretical 
problem  for  the  Mott  transition:  One  must  predict  the  densi 
ties  at  which  excitons  become  unbound  due  to  screening.  A 
theory  for  doing  this  will  be  discussed  in  the  following 
chapter  and  will  be  extended  to  take  conduction  band  anisot¬ 
ropy  into  account. 


III.  Theorv 


As  was  mentioned  in  the  introduction,  Norris  and  Bajaj 
(Norris,  1982)  have  developed  a  theory  for  the  exciton- 
plasma  Mott  transition  in  Si  which  is  in  good  agreement 
with  experiment.  They  obtained  the  exciton  Hamiltonian  by 
assuming  isotropic  masses  and  by  assuming  static  electron-* 
hole  (E-H)  screening  in  the  random  phase  approximation 
(RPA) .  They  then  associated  the  Mott  transition  at  a  given 
temperature  with  the  E-H  pair  density  for  which  the  binding 
energy  of  the  exciton  becomes  zero.  The  binding  energy  was 
evaluated  variationally.  In  this  chapter,  the  above  theory 
will  be  presented  and  extended  to  take  into  account  the  con 
duction  band  anisotropy.  As  was  mentioned  earlier,  this  is 
being  done  in  order  to  assess  the  effects  of  the  electron 
anisotropy  on  the  exciton-plasma  Mott  transition  both  in 
Si  and  (more  importantly)  in  Ge. 


Background 

Theoretical  work  on  the  exciton-plasma  Mott  transition 
was  preceded  by  investigations  of  the  simileu:  problem  of  an 
electron  bound  to  a  donor  impurity  in  a  many-valley  semi¬ 
conductor.  In  the  latter  case,  the  task  is  to  compute  the 
donor  concentration,  N  ,  at  which  the  electrons  become 
unbound  due  to  screening.  Since  the  theory  to  be  presented 
in  this  work  is  a  direct  product  of  the  earlier  work,  the 
donor  impurity  problem  will  now  be  reviewed. 
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Initially,  the  screening  was  taken  into  account  by 
assuming  the  electron  to  be  bound  in  a  Yukawa  potential. 


and  the  variational  calculations  were  done  using  hydrogenic 
trial  functions  (Mott,  1949).  Later,  Rogers  et  al  (Rogers, 


1970)  numerically  integrated  Schrodinger ' s  equation  for 

« 

the  above  case.  The  two  results  were  not  in  very  good 
agreement,  and  hydrogenic  wave  functions  were  seen  to  be 
poor  trial  functions. 

Later,  Lam  and  Varshni  (Lam,  1971)  did  the  variational 
calculation  for  the  Yukawa  potential  using  eigenfunctions 
of  the  Hulthdn  potential.  Their  results  were  in  good  agree¬ 
ment  with  the  calculation  of  Rogers.  The  Hulthdn  potential 
and  its  s-state  eigenfunction,  as  given  by  Greene  et  al 
(Greene,  1977)  are 


2  -vr 

V„(r)  =  -  ^ 


eod  -  e 


) 


(1) 


and 


♦  (£)  = 


1 

2(7ra)^ 


^-(1  -  y/2)r/a  _  ^-(1  +  y/2)r/a 
r 

(2) 


where  e  is  the  electronic  charge,  is  the  static  dielectric 
constant,  a  is  the  first  Bohr  radius  of  the  electron,  and  y 
is  taken  as  a  variational  parameter.  Greene  et  al  used  the 
Hulthdn  wave  function  to  solve  the  problem  of  an  electron 
in  the  Lindhard  potential  (RPA)  at  T  =  0  and  of  an  electron 
in  the  Hubbard-Sham  potential  (which  includes  first  order 


corrections  to  the  RPA  at  T  =  0) .  Their  result  for  the 
Rubbard-Sham  potential  was  in  good  agreement  with  the  result 
of  Martino  et  al  (Martino,  1973) ,  who  solved  the  corresponding 
Schrodinger  equation  numerically.  Thus  the  Hulthdn  wave 
function  was  again  seen  to  be  a  good  trial  function. 

In  all  of  the  above  Vl?ork,  the  electron  masses  were 
taken  to  be  isotropic.  Since  the  electron  masses  in  many- 
valley  semiconductors  are  anisotropic,  Aldrich  (Aldrich,  1977) 
treated  the  problem  of  an  electron  with  anisotropic  mass 
bound  in  the  Lindhard  and  Hubbard-Sham  potentials  at  T  =  0. 

He  performed  the  variational  calculation  using  a  modified 
form  of  the  Hulthdn  wave  function,  which  will  be  used  in 
the  present  work: 


where  a  is  the  Bohr  radius,  S  and  y  are  variational  para¬ 
meters,  and  where  P  is  given  by: 


in»  2 
a^x  + 


5  ^ 


(4) 


Ellipsoidal  energy  surfaces  are  assumed  so  that  longitudinal 
and  transverse  masses,  m^  and  m^,  may  be  introduced,  and 
m*  =  (m^m^)  .  The  parameters  a^  and  a^.  are  given  by 


(5a) 
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where  e  is  a  third  variational  parameter.  The  variational 
parameter,  y,  reflects  the  strength  of  the  Hulth4n  potential, 

B  effectively  expands  or  contracts  the  wave  function,  and  e 
adjusts  the  trial  function  anisotropy. 

The  above  discussion  shows  how  the  choice  of  RPA 
screening  has  come  about  and  how  the  variational  approach 
has  developed.  The  theory  of  Norris  and  Bajaj  for  the 
exciton-plasma  Mott  transition  will  now  be  presented  in 
detail  via  its  extension  to  include  electron  anisotropy. 

This  will  be  done  by  considering  first  the  exciton  Hamil¬ 
tonian,  and  then  the  variational  calculation  for  the  excitonic 
binding  energy. 

The  Exciton  Hamiltonian 

The  exciton  Hamiltonian  is  tedcen  to  be 


+  U(r^ 


(6) 


where  and  (i  =  1,2,3)  are  the  diagonal  elements 

of  the  electron  and  hole  effective  mass  tensors.  The  posi¬ 
tion  vectors  for  electrons  and  holes  are  ^  and  r.  ,  while 
*ei  ^hi  components  of  these  vectors.  U(r^-r^^) 

is  the  interaction  potential. 

In  writing  the  kinetic  energy  term  for  the  electron 


in  equation  (6) ,  ellipsoidal  energy  surfaces  have  been 


asstuned.  The  hole  mass  will  be  taken  to  be  isotropic:  The 
effective  mass  tensor  notation  has  been  retained  in  the 
kinetic  energy  term  for  holes  in  order  to  make  the  deriva¬ 
tions  easier.  The  assumption  of  isotropic  hole  masses  has 
been  made  because  any  attempt  to  go  beyond  this  assvunption 
would  make  the  problem  intractable.  For  a  complete  treat¬ 
ment  of  the  hole  kinetic  energy  see  (Lipari,  1971) ,  and  for 
a  discussion  of  semiconductor  band  structure  see  (Rice, 

1977 :  page  5) . 

A  center  of  mass  transformation  is  now  made  in  order 
to  treat  the  e-h  pair  mathematically  as  a  single  particle 
in  the  asstimed  potential.  Thus 


(7a) 


n» 


ei^ei  ”hi^hi 
%i  +  iHhi 


(7b) 


(7c) 


where  the  m^  are  reduced  masses.  Ignoring  the  translational 
term,  which  does  not  affect  the  binding  energy,  one  obtains 


H 


3 

£ 

i=l 


i  +  U(r) 

"*i  ^ 


(8) 


Since  the  energy  surfaces  are  assumed  ellipsoidal, 
longitudinal  and  transverse  effective  masses  may  be 
introduced : 
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a  term  which  represents  the  average  response  of  the  screen¬ 
ing  carriers  to  the  total  potential.  For  electron  screening 
this  is 

V(r)  =  V^(r)  -  1  - - - —  d^r'  (12) 

y  |r-r'| 

« 

where  V(r)  is  the  total  potential  at  r,  V^(r)  is  the 
unscreened  potential  at  r,  and  <n^^>  is  the  expected 
value  of  the  screening  particle  density  at  r"*. 

The  particle  density  <n^^>  is  given  by 

<n(r)>  =  Trace (p [V]n (r) )  (13) 

where  p[V]  is  the  density  matrix  and  n  is  the  particle 
density  operator.  The  density  matrix  is  obtained  from  its 
equation  of  motion 

-ih||=lS,H]  (14) 

idiere  H  is  the  system  Hamiltonian.  The  density  matrix  and 
Hamiltonian  are  written  in  terms  of  perturbations: 


»  5^  +  «P 

(15a) 

=  H  +  V 
o 

(15b) 

The  subscript  "zero"  refers  to  the  problem  in  the  absence 
of  an  external  potential  and  V  includes  both  the  external 
and  screening  potentials.  The  RPA  arises  when  equation  (12) 
is  assumed  and  the  term  in  5pV  is  dropped  from  equation  (14) 


When  the  £Ubove  analysis  is  carried  out  for  both 
electrons  and  holes,  the  potential  energy  is  foxind  to  be 
given  by 

’'‘a’  ■ '  f?  ^ 

* 

where  e  is  the  electronic  charge,  is  the  static  dielectric 
constant,  ^  is  the  wave-vector  (of  magnitude  g) ,  and  where 
e (£)  is  the  dielectric  function.  The  dielectric  function 
is  given  by 

2 

e(g)  =  1  -  g  (g;)  +  v.g^(2)  1  (17) 

where  and  are  the  band  degeneracy  factors  for  electrons 
and  holes,  and  where  the  g  functions  are  the  density-density 
response  functions  for  electrons  and  holes.  The  response 
function  for  electrons  is 


Mow,  the  integral  in  equation  (18)  has  been  evaluated 


for  practical  computation  by  Meyer  (Meyer,  unpublished) , 
assuming  spherical  energy  surfaces.  In  order  to  do  the 
present  work,  it  was  therefore  necessary  to  assume  ellip¬ 
soidal  energy  surfaces  according  to 


+ 


(19) 


and  cast  the  expression  for  the  dielectric  function  into  a 
fom  where  Meyer's  results  could  be  used.  The  details  of 
this  process  are  outlined  in  the  appendix.  The  result  is 


c(^)  =  1  + 


4irne 

®o^b 


<3(Xgo,ng)  + 


F.  (nv) 


G(X, 


ho 


(20) 


where  n  is  the  electron-hole  pair  density,  the  n's  are  the 
reduced  Fermi  energies  (the  chemical  potentials  divided  by 
kj^T) ,  and  where 


G(x^,ti) 


(irxo) 


^  F=¥TTnr 


I 


f  (x-n)  -^n 


1  + 


1  - 


dx  (21) 


F|,(n)  =  ET  /"x^f(x  -  n)dx  (Fermi  integral;  order  k)  (22) 


m*k.T  L”e^  ^  “et  ^ 


with  a  similar  expression  for  Xj^^.  In  equation  (23) 

®®e  ^  ^®ef”*et^  ^  Fermi-Dirac  function. 

Meyer  has  cast  the  function  G(x^,n)  into  approximate  analy 
tical  form  and  his  results  can  (and  will)  be  used  in  the 
present  work  provided  equation  (23)  is  used  for  x^^  along 
with  the  corresponding  expression  for  Xj^^. 

The  explicit  expression  for  the  potential  energy  term 
in  the  Hamiltonian  has  now  been  obtained  (equation  16  com¬ 
bined  with  equation  20) .  The  variational  calculation  for 
the  exciton  binding  energy  v;ill  now  be  presented. 


The  Variational  Calculation 

It  was  mentioned  at  the  beginning  of  this  chapter,  that 
the  Mott  transition  at  a  given  temperature  is  associated 
with  the  E-H  pair  density  for  which  the  exciton  binding 
energy  is  zero.  The  binding  energy  is  computed  variation- 
ally  by  minimizing  the  expectation  value  of  the  Hamiltonian, 
<H> ,  where  <H>  is  computed  using  the  Hulth4n  wave  function 
given  in  equation  (3) .  This  section  will  therefore  address 
the  determination  of  <H>  «  <T>  +  <V>  ,  where  <T>  and  <V> 
are  the  expectation  values  of  the  kinetic  and  potential 
energies,  and  will  conclude  with  a  discussion  of  the  mini¬ 
mization  process. 


The  Kinetic  Ener 


As  previously  discussed,  the  trial  function  for  the 
variational  calculation  is 


♦<£>  -  ^  (^) 


(24) 


where  r.  is  given  by 


[“/  2  “t 


2  +  *2, 


2  ♦  2 

and  where  a  =  e  4r  /m  e  is  the  exciton  Bohr  radius.  The 

o 

masses  and  are  the  longitudinal  and  transverse  reduced 
masses  introduced  in  the  discussion  of  the  Hamiltonian, 
while  m*  =  (m^m^)^'^^  .  Also,  a^  =  and.  a^.  * 

as  previously  defined.  The  variational  parameters  are  y ,  Br 


and  e. 


The  expectation  value  for  the  kinetic  energy. 


can  be  evaluated  analytically.  The  result  is 


<T>  -  (a^  +  2a^) (4  -  y^) 


where  H*  ■  m*e^/c^4i^  is  the  reduced  mass  Hartree,  the 


energy  unit  used  in  this  work. 


where  V(c[)  is  the  statically  screened  Coulomb  potential 
given  by  equations  (16)  and  (20) . 

It  will  be  seen  later  that  the  determination  of  <U> 
Involves  the  numerical  evaluation  of  a  double  integral  on 
the  unit  square.  It  turns  out  that  if  the  potential,  V(^)  , 
Is  used  "as  is"  in  equation  (29)  ,  then  the  integrand  will 
have  a  finite  value  on  the  side  of  the  unit  square  which 
corresponds  to  infinite  wave-vector  magnitude.  This  diffi¬ 
culty  can  be  eliminated  if  V(^)  is  expressed  as  the  sum  of 
a  "screening"  term  and  a  Coulomb  term.  In  this  case,  the 
expectation  for  the  Coulomb  term  can  be  evaluated  analyti¬ 
cally,  while  the  integration  for  the  screening  term  involves 
an  integrand  which  vanishes  on  the  previously  mentioned  side 
of  the  unit  square. 

Thus  one  writes 
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where 


and 


Vc(2) 


(31) 


(32) 


and  where  £ (^)  is  the  dielectric  function. 

When  equation  (30)  is  substituted  into  equation  (29) 
for  <U>,  two  integrals  result. 

The  Coulomb  integral  is  found  to  be  given  by 

<U^>  =  H*B^g(y)h(c)  (33) 

where  H*  is  the  Hartree,  where  y,  3,  and  e  are  the  varia¬ 
tional  parameters  in  the  trial  wave  function ,  emd  where  the 
functions  g(y)  euid  h(c)  are  given  by 


For  the  screening  term,  one  has 


em^ 


1 


1 

(34b) 


<U  > 
scr 


) 


d^r 


(35) 


Now,  can  be  expressed  as  a  single  triple  integration 

in  q-space  by  making  the  change  of  variables  from  3.  to 
and  interchanging  the  order  of  integration.  One  then  obtains 
the  integral  over  all  q-space  of  (g)  times  the  Fourier 
transform  of  the  wave  function  squared: 


<U  > 
scr 


1 


jF(4»2)Vg  (g)d^q 


(36) 


The  Fourier  transform,  F,  is  easy  to  compute  and  is 
given  by 


where 


(37b) 


-  r  m*  2  .  m*  ,2  .  2  J  ** 

The  masses  which  appear  in  equation  (37b)  are  the  reduced 
masses,  a  is  the  Bohr  radius,  emd  the  quantities  a^  and 
a^  were  defined  in  the  section  where  the  wave  function  was 
introduced  (see  equations  3-5) . 

The  integral  in  equation  (36)  is  evaluated  after  several 
changes  of  variables. 

First,  the  spherical  coordinates  (q,  0,  (j>)  defined  below 
are  introduced. 


=  qcos0 

(38a) 

=  qsin9cos(^ 

(38b) 

=  qsin0sin({> 

(38c) 

The  integration  with  respect  to  can  be  done  immediately, 
and  one  is  left  with  a  double  integral.  However,  because  of 
the  complicated  nature  of  the  integrand,  neither  the  q-  nor 
the  0 -integration  can  be  perfojnned  analytically  and  the 
double  integral  must  be  computed  numerically. 

The  second  and  third  changes  of  variables  in  equation 
(36)  are  to  introduce  respectively  q  =  aq  and  u  =  cos0  . 
Again,  a  is  the  Bohr  radius,  and  the  "tilde"  is  used  on  q 
to  indicate  that  q  is  a  dimensionless  wave-vector  magnitude. 


Thus,  if  the  expression  for  Vg(2^)  in  terms  of  the  di¬ 
electric  function  (see  equation  (31)  is  substituted  into 
equation  (36) ,  one  obtains 


=  I  H*/^du/“dq  jF3(u,q) 


e3(u,q)  -  1 


e3(u 


,^)  j 


The  explicit  expressions  for  F3  and  e3  are  complicated  and 
will  be  presented  after  one  final  change  of  variables.  Also 
subscripts  are  used  on  F  and  e  to  denote  the  slightly  dif¬ 
ferent  functions  which  come  about  because  of  the  changes  of 
variables.  The  subscript,  3,  in  equation  (39)  indicates  the 
third  change  of  variables. 

The  final  change  of  variables  is 


which  tremsforms  the  integral  with  respect  to  q  into  an 
integration  oh  the  interval  0  1.  s  1  .  One  now  obtains 


<U  > 
scr 


2  X  1  1  (  ^c  e4(u,s)  -  1  I 

=  j.  H  /  au/  ds  (41) 


The  explicit  forms  for  F3  and  C3  are  given  below,  while 
the  explicit  forms  for  F^  emd  are  obtained  from  those  for 
Fx  and  by  replacing  q  by 


One  has 


,(u,q)  «  +  tan  -2tan"^(^jj 


with 


?  1^)  ["-("-  "'1 


(43a) 


(43b) 


with 


EjCu.q)  =  1  +  4Mna3)^|j,-^  G(Xeo.ne) 


*■-15  ^'^h^ 

*  F,,^h)  °'*hO 


'%>)  ^ 


*eo  "  I 


5ifer[‘-i‘-fe)-’j’^ 


(44a) 


(44b) 


_  1  m*  H*  l^'hz]  .  “ht]  2  . 


(44c) 


As  introduced  in  equation  (22) ,  the  fractionally-sub- 
scripted  F  functions  are  Fermi  integrals,  and  the  G  function 
is  the  integral  which  is  evaluated  by  Meyer's  interpolation. 

It  should  be  noted  that  all  quantities  in  equations 
(43)  and  (44)  have  been  put  into  forms  where  units  cancel. 
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The  quantity  (na^)  may  be  identified  as  a  dimensionless 

* 

electron-hole  pair  density,  H  may  be  defined  as  a 

dimensionless  reciprocal  temperature,  and  the  quantities 


5i-#^and 

m* 

e 


m*  H* 

lC|_T 
m  D 
h 


may  be  defined  as  dimensionless  electron  and  hole  reciprocal 
temperatures,  respectively. 

Minimization  of  the  Total  Energy 

The  routine  used  to  find  the  binding  energy,  <H> , 
(Shankland,  private  communication)  minimizes  a  function  of 
a  given  number  of  parameters.  In  this  case,  <H>  is  a  func¬ 
tion  of  y ,  $ ,  and  e ,  the  parameters  in  the  trial  wave  func¬ 
tion  introduced  in  Chapter  III.  The  search  for  the  minimum 
involves  random,  gradient,  average,  and  jump  steps,  which 
are  made  initially  with  user  defined  frequencies,  and  finally, 
with  frequencies  generated  in  the  course  of  the  calculation. 

As  convergence  is  obtained,  the  step  size  is  automatically 
reduced  to  insure  that  the  minimum  is  found.  Finally,  jump 
steps  are  made  to  allow  for  the  possibility  of  local  minima. 

There  are  two  possible  approaches  which  can  be  taken  to 
find  the  Mott  transition  density,  is  to  keep  the 

electron-hole  pair  density  as  an  input  parameter  and  minimize 
<H>  as  a  function  of  p,  6,  and  c.  The  density  is  then 
adjusted  manually  until  that  density  is  found  for  which 
^**^minimized  “  ®  other  approach  is  to  include  the 
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density  as  a  variational  parameter  in  the  minimization 

2 

routine  and  minimize  <H>  .  The  first  approach  has  been 
adopted  here  because  it  is  then  easier  to  monitor  the  course 
of  the  calculation. 

The  procedure  is  to  make  a  series  of  short  calls  (25 
function  evaluations  each')  to  the  minimization  routine, 
looking  for  negative  values  of  <H>.  Since  the  goal  is  to 

find  the  density  for  which  "'H^minimized  zero,  and  since 
the  Mott  transition  is  due  to  screening,  one  knows,  as  soon 
as  a  negative  energy  is  obtained,  that  the  density  must  be 
increased.  The  short  calls,  however,  do  not  allow  for  the 
search  step  size  to  be  reduced  below  about  10~^  so  that  this 
initial  search  is  a  rough  one.  Once  a  series  of  positive 
energies  is  obtained,  than  a  long  call  (250  function  evalua¬ 
tions)  is  made  to  check  for  the  best  minimum,  in  this  way 

-4 

the  step  size  is  reduced  to  about  10 

Results  will  now  be  presented  for  Si  and  Ge. 


IV.  Results 


The  material  parameters  used  in  this  work  are  given  in 

Table  I  (See  Rice,  1977;  page  8).  Masses  are  given  in  units 

of  the  electron  free-space  mass.  The  longitudinal  euid 

transverse  electron  masses  are  given  for  Si  and  Ge.  The 

% 

electron  masses  used  in  the  isotropic  approximation  are 
obtained  from  the  optical  mass  average 


m 


eo 


(45) 


while  the  hole  masses,  mj^^  and  reciprocal  of 

the  Dresselhaus-Kip-Kittel  A  parameter  given  by  Rice.  The 
density  of  states  masses,  and  n^h^os'  used  in  the 

calculation  for  the  reduced  Fermi  energies  and  are  deter¬ 
mined  from  the  following  equations: 


-  /  2.1/3 

“edos  =  ^”*e£"‘et^ 


®^dos 


3/2  .  „3/2  \  2/3 
+  m^T  I 


hL 


(46a) 


(46b) 


The  equation  for  the  hole  density  of  states  mass  involves 
the  heavy  and  light  hole  masses,  and  and  thus 

takes  into  account  the  fact  that  the  heavy  and  light  hole 
bands  have  been  replaced  by  a  single  doubly -degenerate  band. 
Table  I  also  gives  the  static  dielectric  constant  (e^)  and 
the  degeneracy  factors  (Vg,V}j)  for  Si  and  Ge. 


The  results  of  Table  I  for  isotropic  Si  are  in  good  agree¬ 
ment  with  the  results  of  Norris  and  Bajaj  (Norris,  1982). 
Only  two  points  have  been  determined  for  anisotropic  Si 


Table 

II: 

Mott 

Transition  in  Si 

DENSITY  (xlO 

cm  ) 

Si  - 

ISOTROPIC 

Si  - 

ANISOTROPIC 

TCK) 

Lower  Bound  Upper  Bound 

Lower  Bound  Upper  Bound 

10’^ 

3.7 

3.8 

« 

7.81 

4.0 

4.1 

4.1 

4.3 

15.6 

5.1 

5.2 

23.4 

6.3 

6.4 

30.0 

7.4 

7.5 

7.8 

8.0 

46.9 

10.1 

10.2 

62.5 

12.7 

12.8 

78.1 

15.2 

15.3 

since 

two 

points 

are  sufficient  to 

show  the 

effects  of  con- 

duction  band  anisotropy.  It  is  seen  that  the  Mott  transi¬ 
tion  shifts  to  higher  densities.  The  shift  is  4%  at 
T  =  7.81”K  and  6%  at  T  =  30.0**K.  These  results  will  be 
Interpreted  after  the  results  for  Ge  are  presented. 

The  results  for  Ge  are  shown  in  Table  III.  Again,  the 
Mott  transition  shifts  to  higher  densities  when  conduction 
band  anisotropy  is  taken  into  account.  Here  there  is  no 
•apparent  shift  at  S^K  while  the  shift  at  the  higher  tem¬ 
peratures  is  approximately  7%. 

The  observations  which  are  to  be  made  here  are  as 
follows:  (1)  The  electron  anisotropy  causes  the  Mott 


Table  III,  Mott  Transition  in  Ge 

DENSITY  (xlO^^cm"^) 

Ge  -  ISOTROPIC  Ge  -  ANISOTROPIC 


T(*K) 

Lower  Bound 

Upper  Bound 

Lower  Bound 

Upper  Bound 

5.0 

.13 

.15 

.13 

.15 

12.5 

.26 

• 

.27 

.28 

.29 

20.0 

.39 

.395 

.41 

.43 

tramsition  to  shift  to  higher  densities.  (2)  The  shift  is 
greater  for  high  temperatures  than  for  low  temperatures. 

(3)  The  shift  appears  to  be  no  greater  in  Ge  than  in  Si. 

The  first  observation  can  be  explained  by  the  fact  that 
the  anisotropy  will  reduce  the  ability  of  the  electrons  to 
screen  emd  thus  a  higher  density  will  be  required  to  screen 
out  the  Coulomb  potential.  The  second  observation  can  be 
partially  explained  by  noting  that  the  electrons  and  holes 
have  more  thermal  energy  at  higher  temperatures  and  energetic 
carriers  do  not  screen  as  effectively  as  less  energetic 
carriers.  The  third  observation  is  partially  explained  by 
a  calculation  of  the  exciton  binding  energy  in  the  cdssence 
of  screening. 

Table  IV  shows  binding  energies  in  Ge  for  an  electron 
bound  to  a  hole  and  to  a  donor  impurity r  in  the  absence  of 
screening.  It  can  be  seen  that  the  anisotropy  does  not 
significantly  affect  the  exciton  binding  energy.  This  is 


Table  IV;  Binding  Energies  for  Zero  Screening  (Ge) 


Exciton 

Isotropic  Masses 
Anisotropic  Masses 

Donor  Impurity 

Isotropic  Masses 
Anisotropic  Masses 


Binding  Energy  (meV) 
2.67 
2.73 

6.91 

9.78 


because  the  reduced  masses  are  dominated  by  the  light  masses 
and  in  Ge  one  has  m^^  =  .082  and  mj^  =  nij^^  *  .075  .  Thus 
the  ratio  of  longitudinal  to  transverse  reduced  mass  is  on 
the  order  of  only  two  for  the  exciton.  For  the  donor 
impurity,  the  above  situation  no  longer  prevails  and  the 
anisotropy  significantly  affects  the  binding  energy.  The 
result  of  9.78  meV  is  in  good  agreement  with  Faulkner 
(Faulkner,  1969)  who  obtained  9.81  meV  for  an  electron 
bound  to  a  donor  impurity  in  Ge. 

Thus,  the  shift  of  the  Mott  transition  densities  is 
due  mainly  to  the  effects  of  the  anisotropy  on  the  screening 
since  the  anisotropy  does  not  significantly  affect  the 
exciton  binding  energy. 


V,  Summary  and  Conclusions 

In  this  work,  conduction  band  anisotropy  has  been 
incorporated  into  the  theory  of  the  exci ton-plasma  Mott 
transition  (Norris,  1982) .  The  exciton  was  treated  mathe¬ 
matically  as  a  particle  in  a  screened  Coulomb  potential 

« 

where  static  electron-hole  screening  in  the  randon  phase 
approximation  was  assumed.  The  Mott  transition  was  asso¬ 
ciated  with  the  electron-hole  pair  density  at  which  the 
exciton  binding  energy  in  the  assumed  potential  is  zero, 
and  the  binding  energy  was  computed  variationally  using  the 
ground  state  eigenfunction  of  the  Hulthdn  potential  as  a 
trial  function. 

The  results  obtained  from  the  above  theory  lead  to  the 
following  conclusions: 

(1)  The  conduction  band  anisotropy  does  not  signifi¬ 
cantly  affect  the  exciton  binding  energy. 

(2)  The  conduction  band  anisotropy  decreases  the 
ability  of  the  electrons  to  screen  and  thus  increases  the 
Mott  transition  densities  beyond  those  predicted  by  RPA 
screening  with  the  conduction  band  taken  as  isotropic. 

(3)  The  effects  of  the  electron  anisotropy  are  no 


more  pronounced  in  Ge  than  in  Si 
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endix:  The  Anisotropic  Dielectric  Function 


The  screened  Coulomb  potential  in  wave-vector  space  is 


given  by 


V(a)  =  - 


c.q  €  (gr) 


where  e  is  the  electronic  charge,  is  the  static  dielectric 
constant,  q  is  the  wave-vector  magnitude,  and  e  (g^)  is  the 
Lindhard  dielectric  function. 

The  dielectric  function  is  given  by 


e(a)  =  1  - 


where  and  Vj^  are  the  degeneracy  factors  and  where 
and  are  the  density-density  response  functions  for 

electrons  and  holes  respectively.  The  response  function 
for  electrons  is  given  by 

f  [  (Eg  (k)  -Ug)  AwT]  -f  [E  (k-q)  -y g)  AjjT  ] 


In  equation  (18) ,  y  is  the  chemical  potential  for  electrons 
and  f(x)  =  l/(e*  +  1)  ,  the  Fermi -Dirac  function.  The 
function,  Eg(k},  is  defined  by 


E.  (k)  =  k?  +  (k?.  +  k?) 


where  and  are  the  longitudinal  and  transverse  electron 
effective  masses. 

The  explicit  expression  for  the  dielectric  function  is 
obtained  by  performing  the  integration  in  equation  (18) ,  with 
an  almost  identical  calculation  for  holes.  The  first  step  is 
to  make  the  following  changes  of  variables  in  order  to  make 

^  “9-^  spherically  symmetric  (Aldrich,  1977:  2724): 


■Cf/ 


.  /nisi 


l'^\ 

\<l 


■  rtf 

■m 

■itf 


(A-la) 


(A-lb) 


(A-lc) 


In  equation  (A-1) ,  m*  is  given  by  m*  =  ^™e-£"'et^  * 

With  the  ^d^ove  changes  of  variables,  it  is  found  that 


E_(k)  »  =  E_(s) 

®  “  2m!  ® 


(A-2a) 


2 

(k  -  a)  =  A^|s  -  ti^  =  E^(s  -  t) 


(A-2b) 
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Of  course,  the  E^'s  on  the  two  sides  of  equation  (A-2) 

represent  different  functions,  but  the  arguments  will 

always  be  explicitly  stated,  and  no  confusion  should  result. 

3  3 

It  is  found  that  d  s  =  d  k  so  that  becomes 


fI(Eg(s)-iJ^)/kjjTl 

Ee(s)-Ee(s-t) 


f[  (Eg(S“t)  -yg)/kj^T] 
Ee(s)-Eg(s-t) 


(A-3) 


Both  of  the  integrals  in  equation  (A-3)  can  be  done  in 
polar  coordinates.  For  both  integrals,  the  s_-axis  is  chosen 
to  lie  along  the  vector  t.  Thus,  for  the  first  integral,  one 
has 

2 

E  (s  -  t)  =  -^(s^  +  t^  -  2stcos0)  (A-4) 


For  the  second  integral,  one  makes  the  substitution 
s*  =  s  -  t  ,  whereby  t  is  eliminated  from  the  numerator. 
Since  d  s  =  d  s  ,  the  primes  may  be  dropped  and  one  then 
obtains 

2 

E.(s  +  t)  =  Ar(s^  +  t^  +  2stcose)  (A-5) 


in  the  denominator  of  the  second  integral. 

After  performing  the  anguleu:  integrations  in  equation 


<A-3) ,  it  is  found  that 


V 


(A-6) 


Equation  (A-6)  is  in  a  form  which  begins  to  resemble 

the  result  obtained  by  Meyer  (Meyer,  unpublished) .  The  final 

« 

expression  for  the  dielectric  function  will  be  the  same  as 
Meyer's,  except  that  scaled  arguments  will  be  used. 

The  reduced  Fermi  energy,  introduced 

here.  One  defines 

2  2 

and  makes  the  change  of  variables  x  =  (iS  /2m*kjjT)  s  .  Then 
equation  (A-6)  becomes 


(A-8) 


Now,  equation  (A-8)  does  not  contain  the  electron-hole 
pair  density,  n.  To  introduce  the  density,  it  is  necessary 
to  evaluate 

n  =  fl(Eg(k)  -  Ug)AbT]  (A-9) 

"When  the  integration  in  equation  (A-9)  is  performed,  and  the 
result  is  substituted  into  equation  (A-8) ,  equations  (20) 
through  (23)  for  the  dielectric  function  are  obtained. 
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